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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO. 960 

FORCE AND MOMENT COEFFICIENTS FOR A THIN AIRFOIL 
WITH FLAP AND TAR IN A FORM USEFUL FOR 
STABILITY AND CONTROL CALCULATIONS 

i 

By Roland J. White and Dean 5. Klampe 

SUMMARY 


Recent airplane design trends have teen directed along 
the line of cont r ol- f r e e- st ahi li t y analysis involving motions 
of the control surfaces having frequencies loss than those 
considered in flutter calculations, hut great enough to re- 
quire a knowledge of the air forces resulting from velocity 
displacements. In addition to this, control arrangements 
have been considered having the entire airfoil pivoted or 
having a flap or tab actuated by pressures developed on the 
contour of the airfoil. In order to facilitate these calcu- 
lations a systematic presentation of the necessary pressure 
and force coefficients acting on a thin airfoil having a flap 
and tab has been made and presented in this report in a form_ 
suitable for stability calculations. 

Equations for the pressures have been derived from the 
basic equations of reference 1 and their development has been 
given in detail. Force coefficients have been taken direotly 
from references 1 and 2 or integrated graphically from the' 
curves of pressure distribution. In order to preserve a con- 
tinuity of equations, the displacement, velocity, and accelt- 
eration stability coefficients all have been calculated even 
though the acceleration or virtual mass effect generally may 
be neglectod in stability calculations. In all cases only 
the real part.of the complex equations of reference 1 is con- 
sidered due to the order of frequencies generally involved - . 


RESTRI OTED 
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SYMBOLS 

a angle of attack referred to. relative wind (radians) 

6, a 1 angle of rotation of wing referred to horizontal (radians) 

• — , Y angle between the relative wind vector and horizontal 
v (radians) 

6,$ flap deflection (radians) 

6t tab deflection (radians) 

b semi-chord of airfoil (ft) 

c flap hinge position in fractions of the airfoil semi- 
chord (measured from airfoil center) 

a center of rotation of airfoil in fractions of the airfoil 
semi-chord (measured from airfoil center) 

* 

H center of rotation of airfoil in fractions of the total 
chord (measured from L.E.) 

E ratio of the flap chord to the total chord 

cqi total chord of airfoil' (ft) 

Cf flap chord (ft) 

Cfc tab chord (ft) 

h vertical displacement of airfoil (ft) 

v velocity (fps) 

# 

M.A.O. mean aerodynamic chord (ft) 

x distance from airfoil chord midpoint to any chordwise 
station measured in fractions of the semi-chord 

p pressure difference between upper and lower surfaces on 

airfoil 

t time (sec) 
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p mass air density ( slugp/ft 3 ) 

cp velocity potential function 

r circulation 

P = p^/i pv s 


Cl 

Cm 

C H 


°I () 


lift coefficient 

pitching moment* coefficient 

hinge moment coefficient 
QCj, 

3 () 


C D () = ^ 
U 3 () 

r a °H 

° H o c m 


() = 


a p 

■ ao 


HOT AT I OH 


Before the stability coefficients may be derived from 
the equations of reference 1 it is necessary to distinguish 
between the degrees of freedom employed in reference 1 and 
•those generally used in control-free-stability calculations. 


Oonsidor an airfoil to have its center of motion position 
at point a as shown in figure 1. The angle 0 as shown in 
this figure Trill define the angle of the 'chord with respect to 
the horizontal and the angle V will be the angle of the 
flight path 'with respect to the horizontal. Then 


0 = 0 ,' 



where a* denotes the value of a- as used in reference 1 
to distinguish from the value of a, employed in stability 
calculations. 
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The angle of attack of the airfoil is then 

A , , fi. 

a = 6- 'v = a'+ — 

' v 

• 

where h is the translational velocity defined in reference 1. 
In the stability notation chosen the differential operator 
D() will indicate 3 ( ) /3t where T is the aerodynamic time 
having units of half-chord lengths. The value of t is de- 
fined by the equation Tb = vt , giving 

3t /3 t = v/b 


where 

v true airspeed (ft/sec) 

t time ( sec ) 

b half chord of airfoil (ft) 

The degrees of freedom of the airfoil will be a, 6, 8, and 
8 t where 6 is the flap angle and is equivalent to (3 used - 

in reference 1, and 6 ^ is the tab angle. All angles will 
be expressed in radian measure. 

In forming the stability coefficients each motion is 
assumed to take place with the other three degrees of free- 
doms maintaining constant values. It is now necessary to re- 
place various terms of the equations given in reference 1 with 
the newly defined coordinates. The substitutions are given in 
table I and are all apparent with the exception of those for 
D(8) and D 3 (6). 


The value of D(0) expresses tho airfoil rotating at 
uniform angular velocity but maintaining a constant angle of 
attack a. To accomplish this the flight path angle must have 
a constant rate of change producing a ..curved flight path; 

•hence 


giving 


9 - H-lv - D(0) 


D(6) . ii . is' 

V V 

D s (e) = ^1’ 


V 

¥ 


v 
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If a is to "be constant then 


giving 


hence 


a , 1 = e 


— + constant 
v 


•• 



D(6) £ 


h 

v 


-D(e) 



In order to obtain equations in terms of D(0), a substitu- 
tion for &' and £ from the equations of reference 1 must 
be made . 

In case the airfoil considered is a tail plane the aero- 
dynamic time T will be based upon the main wing chord which 
may be the wing mean aerodynamic chord. In this case the 
valuo of a is to be measurod in terms of half-chord lengths 
of tho tail chord and the final valuo of the stability’ or 
pressure coefficient converted by multiplying by C^/M.A.G, 

for D ( ) and (Cgi/M.A.C.) 2 for D S () coefficients, where 
Orjt is the tail chord. If l is the distance of the aero- 
dynamic center of the tail plane aft of the airplane center 
of gravity in feet, then 


a = 


2 ± _ 1 . 

On? ** 2 


The stability and pressure coefficients have the subscript 
denoting the type of. coefficient;' for example, 


P B (0 ) 


°Hd 3 ( 6 


SP 

" dD (6 ) 

5 C E 

> ~ 3D S (5) 
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TABLE FOR CONVERTING EQUATIONS OF REFERENCE 1 TO FORM 
EQUATIONS OF THE STABILITY COEFFICIENTS 


To obtain an 
equation in: 

Replace the 
following terms 
of reference 1* 

By 

a 

h/v 

a 

D (a) 

bii/ v 2 

D (a) 

3>( G) 

• 

ha 1 /v 
. h£/v s 

D( 0) 
-D(6) 

D 3 (6> 

h 2 a'/v 2 

D 2 (6) 

6 

P 

6 

D(8) 

bp/v 

D(8) 

D 3 (8) 

b 3 £/v 3 

D 3 (fi) 


*A11 other terms of equations taken as zero 

C = 1 ~ E - H 

2 


PRESSURE DISTRIBUTION OVER A THIN AIRFOIL WITH FLAP 


The chordwise distribution of pressure difference' over a 
thin airfoil with flap will he calculated using the equations 
of the non— circulat ory and circulatory velocity potentials 
given in reference 1. It will he assumed. that the frequency 
of oscillation is low enough for terms containing G, as de- 
fined in reference 1, to he taken equal to zero. The part of _ 
the pr e s sur e c di st r ihut i on not containing 2ttF will he defined 
as the n hasi^ pressure di st r ihut i on , ,f and that part containing 
2ttF will define the additional pressure distribution. This 
segregation will permit airfoil characteristics to he calcu- 
lated for airfoils having a finite aspect ratio. From page 6 
of reference 1, the pressure difference p between the upper 
and the lower surface of the airfoil can he written as: 
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A positive pressure difference is taken as upward, which, 
is the reverse of reference 1. From page 5 of reference l t 
the velocity potential of the non- cir culat ory flow, using the 
notation of reference 1 and substituting a 1 for a, is 



( 4 ) 
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wher e 


3 16 g g _ 1 3| 
3-x H 3 s 



If equations (4) and (6) are substituted in equation (l), 
the pressure distribution for the non-circulatory flow can be 
found. It is now necessary to determine such a circulatory 
flow that the pressure between upper and lower surfaces at the 
trailing edge is zero. If this pressure distribution is added 
to the above distribution, then the pressure distribution duo 
to the movement of the airfoil having an infinite aspect ratio 
can be found. Pressures determined by this method will pro- 
duce a lift; hence the angle of attack of the airfoil will be 
reduced so that the net lift is zero by deducting the pressure 
distribution due to a 1 in proportion to the effective change 
in angle of attack. The resulting pressure '"di stribution then 
will be termed the "basic pressure distribution" and will be 
independent of aspect ratio. 

The pressure distribution due to the non-circulatory flow 
obtained by substituting equations (4) and (6) in equation (l) 
is 
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Po = - V s 
2 


-4 s: 

^ J 1 — 2 


a' + £- v 
2 


2 \ 


2(3 +' 2ax - 4x 3 ) 


ba, 1 


tj 1 ** 


£ T =[ 4 d-a)/rr^]^ ♦ 


£ v a 


~4x 


L yrr 


h 

v 


+ £ v 3 
2 


4^1 ~ 




M'+ £ v 3 

V s 2 


4_ ^ -x cos" 1 c 

TT 


yrr 


- log N 


(x - c) 5 log IT 
dx 


6 + £ V s 

2 



*1 -x A 
v ' 1 - x‘ 


(3 + 2cx - 4x 3 ) _ x k _ „ 

+ cos c - 4(x~ c) log IT 

•J 1 - x 2 


.(x - c) 3 Q..A° 6 * 
dx 


M. + £ v 3 
^ 2 


£( yi - *vi « 


- 2c cos -1 c + x cob"* 1 c) - ( x - ©)* log N 


-1 . 2 


U- (7) 
2 


The .pressure distribution due to the circulatory flow neg- 
lecting the wake effects can be found by letting x 0 — i>c» in 
the equation for (2 tt/AP) dcp/dx as given on p_age 6 of refer 
ence 1* Then,when AP = equation (l) gives 




( 8 > 


Now for zero pressure at the trailing eige 


Po + Pp = 0 


for x 


1 


(9) 
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In order to avoid an indeterminate condition at x = 1, mul- 
tiply p 0 and p,p . "by -/l - x 2 before substituting in equa- 
tion (9 ) . Then 



*• 


X* 


as defined in reference 1, If equation (10) is substituted 
in equation (8), the pressure distribution becomes 

a \ Mi + * + p 

/V V TT 

+ iii ill (12) 

2tt v _ 



The total pressure distribution written in terms of the pres- 
sure ratio P for an airfoil of infinite aspect ratio is 


p = (p o + PrO/§ v z 
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To form the basic pressure distribution, the angle of attack 
a* will be altered so that -T becomes zero. The pressure 
distribution due to a change in a* only is from equation 
(13) 

P « = 4 j | " ■ -£ a' (14) 

a v 1 + x 

Now, whenever a virtual or actual t change in the airfoil 
camber occurs due to a', 6, or 6, the angle a,' of equa- 
tion (10) must be changed so that r i s zero. This must be 
done since by definition the basic pressure distribution is 
selected for zero lift. For this condition the following 
pressure distribution must be deducted from equation (13) to 
farm the basic pressure distribution. 




T 


1 o 


TT 


+ 


Ell 

2t r 



then the basic pressure distribution which will occur for 
r = 0 is 

PB = P - Pr 0 



8 yi-x 3 - 2 /isS. \ 

ba 1 + f 4 /x 

J 1-x 2 

1 .3" , r 

* a • + 4 



L J i+xj 

v ^ \2 

/ 

I v" l 

- 


r* /f4 

rr^i J 1-x 3 


- log H - (x-c) 


6 + 2 


TT 


(3+j x r 4 x S -l C08 -i c 


L ,/T-x ; 


+ ( 1- c ) x / - 4(x-^c) log H r. (x-c)® - — 


N 


1— x* 


3x 


b§_ 

* v 


^/l-c 3 + (x-2c) cos -1 c^ - (x-c) 3 log 


H 


i-i * (16) 

v 3 


For a wing having a finite span the slope of the lift curve m 
will be less than 2 tt , This will cause P a ' from equation (14.) 

to be reduced by the factor J5_. If the lift effectiveness 

2rr 

due to the change in camber is defined as 
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Kg = lio, 

P TT 


P A - f 

TT , 


1 - X 

1 + X 


then the total pressure ratio for a finite wing having a 
slope of lift curve m can he written as 

P = Lpa] a' + fp^ + P A 1 + ["p B( J' - a P B£ ]^l 

LJ L J T L Jv 3 


+ mP A * + £g£ + P B |3 + P A P + [ P Bg + “ K f3 PA 4 - 

" + [pbi? ^pr (is: 


where 


8 Yl^x 3 - 2 / l-.E 

V 1+x 


= 4 V l-x a 


1 x 



- log N - (x-c) 


= 2l [~ L 3+x~4x 3 )_ 00B -i c + (i_c)x / 


l<-x s 


— 4 ( x- c) log S' - (x-c) - — — — 

3 x 


= -| J 1-x 3 i 


(x-c) log S 
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Yalues of the coefficients defined "by equations ( 17 ), 
( 18 ), and (20) have teen calculated and plotted in figures 3 
to 7 . 


LIM, PITCHING MOMENT, AND FLAP HINGE MOMENT POE A 
THIN AIEPOIL WITH A FLAP 


Prom the equations for the pressure distribution deter- 
mined in the provious section, tho airfoil lift and momont 
coefficients can he obtained by integration. These integra- 
tions have been performed in reference (l) for the case of 
the airfoil with a single flap. Hence, it will only be 
necessary to rewrite the equations here using the necessary 
notations to convort to coefficient form. Let 


Or, = 


Cm - 


°H = 


-P 

p b v‘ 


M, 


a_ 


2 p b 3 v 3 
Me 


!BL 


o 3 «, 3 3 

2 p v b E 


whor o 

P lift force given by equation (XVIII) of reference 1 


M 


a 


pitching moment given by equation (XX) of reference 1 


Mp hinge moment given by equation (XIX) of reference 1 
ratio of flap chord tc total chord 


Mp 

hinge 

E 

ratio 

2 rr 

P = m, 

then 


°l - 

met 1 + 


rr + ( — -a ) m 

b.a.L 4. 

— an 

C5 •• 

b a 1 

r 

+ m 

i + 

TT 

L \2 / 

V 

_ __ 

v s 


V 

_ __ 


bS 

■n.3- 


+ 

m 

8 4- -T 4 + m 1 

bg_ + 



TT 

Stt J 

v j 

— — 


bf£ 

^3 


( 21 ) 
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Cm 
w m 


+ 


+ 


+ 


f m( + JL 

L V. 2 4, 

)]“' 

+ [- i(l - a ) + 

m i.| 
2 ' 


[- l(i + 

“)] 

h S a' + a + 1 

V s L \ 3 4 

w 

+ j aiT j hh 

L 2 j v a 

r (T 4+11 

L 2 

2-1 + n 

<t + 0 ^ s + [ 

-l(’ 

’i-Ta - (c-a)T 4 4 1ll\ 

Si + 

—1 

si + fh. + (c-a) 

hJ, 

ifi ( 22 ) 

V 2 

2tt J 

v 1 2 

2 J 

V 3 



All -values of T are a function of a or c only with 
the exception of T 1S and T g . 

* 1 = = - \ ( T 7 + I I ) ' 

From reference 3 the value of T 17 is 

T t , = -21. — + (a - Cl-c s ) 



For the purposes of this report the following coeffi- 
cients defined in terms of the notation of re.ference 1 will 
he adopted; 
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Ka* 


= (l - 
V2 


k B » 5 io/tt 


°L B p = T4 


° l b£ 


= -I. 


J mA 




OHa = - t is/4ttE' 


Kp 

= T 1 1 / 2tt 

E 

= c f /c T 

°L B i' 

1= IT 

° m Ba' 


°^Ba 1 

1 — S.TT 

C niBa' 

■ - f (t + a °) 

° L Bh 

= TT 

°mBh 

= — a 
2 

° L Bp 

= 0 

C m B p 

= ( T 4 + Iio ) / 2 


->-(24) 


r 1 * 
msp 


= - (T i - *8 


L > 


+ ±i±l/2. 


° m BP " - 


Tj_ 

2 


C HBa' “ - I i?/ 2S 

°Hbp 

°HBa' = 

°Hbp 

CH B £ = T X /2E 3 

° h bp 


= -(T 5 - T 4 T 10 )/2ttB 3 

OHbr = 


Values of these 
in figures 2, 3, 

coefficients have been 
and 8 to 13. Then 

calculated and 

plotted 

j 

0 x, = na 1 + j 

[_° L Bi 

, + mK^iJ + £ 

°LBa'_ 

i v 3 " r v, 

bet,' + m h + 

j 2 V 

"c T -1 Jz£ 
L BhJ v 3 

+ ! 

nKp 

P + [cx, B p + mK p] 

bP + 

V 

Tex, -1 ^ 
[_ Cl bbj v 2 

( 25) 
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1 1 

11 

a 

o 

- — r 

C_ . m 
m A 

' <x‘ + [ 


+ °m A ] 4" 

+ [°®b5'] 

b s a ' 
v s 


+ 

_° *A 

l + [‘ 

3m Bh] 

bh + 

V s 

[ Cm BP 

+ mK p 

°* A ] ? 



+ 

_° m B| 4 

• mKp 0. 

„ ( 
m A " 

c-a ) 
2 

° L Bp] 

b£ 

V 




+ 

L °mB6 ' 

. (c-a) 
2 

°Lb£ 

v 3 





(26) 

o 

H 

II 

[m C Ha 

a, 1 + 


1 + m 

K i' °Ha] 4 

l 



+ 


( c-a 
2 

1 “ + 
J V s 

m 0 Ha 

i + [ oh b£ 

1 bh 

J V 3 


+ 

[°h bp ' 

f m Kp 

°Ha] p 

+ [o 

Hbs + 

m Kp C 

ha]4 + [° 

h bp] 

b 2 p (27 
v a 


SUMMARY OS' STABILITY AND PRESSURE COEFFICIENTS 


Following is a list of equations for calculating the 
pressure and stability derivatives. In some cases equations 
involving S t are not written, since the tab effects can be 

calculated from the equations for 6 by use of the proper 
chord ratios . 
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I . Pressure Ooefficients 


= m P A 




P D(<x) = 

p D(e) = 

p D 3 (0) = 
P* = 


■ 3 h 


J VM.A.C./ 

>■ - p >s + ■ Pa ] kzr) 
'***' - a Pb£ ] (irhrT 

PbP + m p aJ 


p 3>(8) - [ e bp * » *8 Pa] (jr^x) 

*»•(.) -M (ioV:) 2 


II . Airfoil Lift Coefficients 


°La = m 
°I»D(a) es 
°Ld( 6 ) = 
°LD S (e) = 


_°LBh] ( M .a!o.) 

OlBa' - °B 3 £ + "**■] (iXr) 

Qij 
M, 


0be »'J Onx) 

°l e - “S'Pj 

°Bd( 6 ) =[°PB 8 ♦ » K «] (ixx) 




(28) 


( 29 ) 
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III, Airfoil Pitching Moment Coefficients 
°m a = |^ m °mA 




C ^D(a) " (k,a!c.) 

°“D(e) = 


°niBa' ~ °*Bh. + m K a 1 °m A 


] (jtt:) 


G mn 2 (e) = 

C m6 = 

°nn(8) = 


Ci 




lmB a' I Vm.a.c 

® m B{3 + m ^ m A 


i 

V (30) 


>m B |3 - ( 1-J5— H ) 0l 3 | + m Kg 0 mA 


C<p 


:) 


°”D S (6) * °”Bg - °L B f(] ( M .1*07) 

IV . Flap Einge Moment Coefficients 


°Hfa " 
°HfDU) = 
°HfD(6) = 
C HfD 3 (e) = 
°Hf 6 = 
°HfD(6)-> = 
°HfD 3 ( 6 ) = 
°Hf8t “ 
°HfD(8t) = 
°HfD 3 (6t) = 


m C Ha J 

° H sh] (m.a%t) 

- 0 =lh. * m °S A ] (jjx) 

- '—r 1 ° L 3i] („.A%.) 

Ch sp + m Kp C Ha 

g h 3 8 + m K -P °Ha ( m .a T cT) 

0kb p_ (m.a^ct) 

°H58t + m E ' St ° HA 

°HB 6 t + m K H °EaJ („ a T 0 :) 


(31) 
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V. Tab Hinge Moment Coefficients 
°H fca = m 0 Ha 

‘ [° KB4 + " X4 H (idx) 

°Hts t " + m Kp °H A 

°3tD(6 t ) = [ Ce BP + °H A J 

08w * ( '* ) 

°H t6 = C HtB p + m Kp C Ha 

°EtD(8) = GE tB0 +m£ ? °HaJ ^ c -) 

=H tD = (6 , -[<*«?] 

The basic tab hinge moment derivatives OEfcsp » ^HtBS > 
and were obtained from a graphical calculation of the 

first moments about the tab hinge axis of the areas under the 
curves of , Pb( 3* and PBg . These curves are presented 

in figures 14, 15, and 16. 

Curves of Ojpgg^. which are shown in figure 17 are taken 
from reference 2, and the curves of Cn B g t and Cjj^g^ , 

figures 18 and 19, were constructed from a graphical calcu- 
lation of t.he first moments about the flap hinge axis of the 
areas under the curves of P B £ and P-gjjj for the deflected 
t ab . 



Curt i ss-Wr i ght Corporation 

Lambert Pield, Saint Louis, Ko., April 11, 1944. 


NASA TN No. 960 


20 


REFERENCES 


1. Thaodor sen, Th.eod.orei General Theory of Aerodynamic 

Instability and the Mechanism of Flutter. NASA Rep. 
No. 496, 1935. 

2. Ames, Milton B., Jr., and Sears, Richard I.: Determina- 

tion of Control-Surface Characteristics from NACA 
Plain-Flap and Tab Data. NACA Rep. No. 721, 1941. 

3. Theodor sen, Theodore, and Garrick, I. E.: Mechanism of 

Elutter - A Theoretical and Experimental Investigation 
of the Elutter Problem. NACA Rep. No. 685, 1940. 


- * 




NACA TN NO. 960 


Fig- 2 



Figure 2.- Airfoil effectiveness . 
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Figure 3.- Additional force coefficients. 
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Figure 4.- Basic pressure distributions for angle of attack changes. 
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Figure 5.- Basic pressure distribution P B 0 UO, 16. 20, £5, 30, 35, and 40 percent chord flaps) . 
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Figure 6.- Sasio pressure distribution Pap {10, 15, 30, 35, 30, 35, and 40 percent chord flaps). 
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a 


Pbp * § [Vl - x s ( V 1 -c 2 + (x - So) cca^c) 
t 1 tx - o) 2 log h] K 1 f- 


- 1.0 -.8 -.6 


Figure 7.- Basic pressure distribution P B p CIO, 15, 80, 85, 30, 35, and 40 percent chord flaps). 
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Figure 10.- Basic pitching moment coefficient* for angle of attack change*. 
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Figure 11.- Basic- pitohing moment coefficients for elevator angle changes 
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“ flap chord 

Figure 14.- Basic tah hinge moment derivative, , 



for elevator changes (per radian). 
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Fig. 14 






NACA TN No. 960 Jig. 15 




I 



_ tab chord 
cj flap chord 

figure 16.- Basic tab hinge moment derivative, Ctt for elevator changes 
(per radian). 
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c f flap chord 
Cjj " airfoil chord 


Figure 17.- Basic elevator hinge moment derivative 
tah changes (per radian) . 














Cf flap chord 
Cg, airfoil chord 


Figure 19.- Basic elevator hinge moment derivative, 
changes (per radian) . 
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